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Abstract 

A coupled model containing two neurons and one astrocyte is constructed by integrating Hodgkin-Huxley neuronal model 
and Li-Rinzel calcium model. Based on this hybrid model, information transmission between neurons is studied numerically. 
Our results show that when the successive spikes are produced in neuron 1 (N1), the bursting-like spikes (BLSs) occur in two 
neurons simultaneously during the spikes being transferred to neuron 2 (N2). The existence of the astrocyte and a higher 
expression level of mGluRs facilitate the occurrence of BLSs, but the rate of occurrence is not sensitive to the parameters. 
Furthermore, time delay % occurs during the information transmission, and t is almost independent of the effect of the 
astrocyte. Additionally, we found that low coupling strength may result in the distortion of the information, and this 
distortion is also proven to be almost independent of the astrocyte. 
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Introduction 

Although the number of the glial cells is several times larger 
than that of the neurons in most parts of the brain, few studies 
have focused on the effect of glial cells on neuronal behavior. Over 
the past decades, an increasing number of works have demon- 
strated that the interaction between glial cell and neuron serves an 
important function in information transmission in the neuron 
system [1-3]. Astrocytes are the most numerous type and the best 
studied glial cells. Astrocytes modulate synaptic transmission 
through many different pathways [4,5]. The most discussed one 
is that the presynaptic neuron release a kind of neurotransmitter, 
glutamate, which activates glutamate ionotropic receptors (i- 
GluRs) on the postsynaptic membrane. Astrocytes participate in 
this synaptic transmission by responding to the glutamate in the 
synaptic cleft through calcium elevation; this elevation of Ca 2+ 
above a certain threshold triggers the release of glutamate to the 
synaptic cleft [6-9]. Calcium elevation also results in the release of 
other transmitters such as Adenosine Triphosphate (ATP). ATP 
activates the purinergic ionotropic receptors, which facilitates the 
enhancement of neuronal excitability [10-15]. Moreover, astro- 
cytes absorb excess potassium released by neurons in the synaptic 
space and thus regulate excitation [16,17]. This bidirectional 
coupling between neurons and astrocytes indicates the concept of 
the "tripartite synapse" [18-21]. 

Traditional modeling studies of neuron system consider the 
coupling between the neurons, but ignore the participation of glial 
cells [22-25]. Nadkarni and Jung introduced a model accounting 
for the interaction between the neurons and the astrocytes [26]. 
They model the effect of astrocyte on neuron through a calcium- 
dependent inward current in the neuron. The calcium-dependent 
function is fitted from experimental data [27]. This kind of 



modeling scheme is extensively employed by many researchers 
[14,15,20,28,29]. The model proposed by Nadkarni and Jung 
predicted the seizure-like spontaneous oscillations in the absence of 
stimuli. Following Nadkarni and Jung, many modeling studies 
focus on the contribution of astrocytes to epilepsy [30-34] 
motivated by experimental findings [35]. For example, Amiri et 
al. concluded that disruption of the homeostatic function of 
astrocytes may initiate the hypersynchronous firing of neurons 
through successive research works [32-34]. This finding suggests 
that the neuron-astrocyte interaction may represent a novel target 
to develop effective therapeutic strategies for epilepsy. 

Neurons are widely accepted to be organized into networks, and 
neuronal networks exchange information through electrical and 
chemical synapses. Increasing evidences indicate that astrocytes 
are also organized into networks [4], and astrocyte networks are 
interconnected through gap junction channels. The channels are 
regulated by extra- and intracellular signals that enable the 
exchange of information. Based on these two networks, a recent 
review paper suggests the concept of "astroglial networks" [2]. 
Many recent modeling works focus on the neuronal synchroniza- 
tion in the astroglial network [28,36 — 40] . As an example, by 
integrating Norris-Lecar neuron model and Li-Rinzel calcium 
model, Amiri et al. constructed a model to study how astrocytes 
participate in the interplay between the pyramidal cells and 
interneurons [40]. Furthermore, they extended their three-unit 
model to a neuronal population model to study the effect of 
astrocyte on neuronal synchronization. Astrocytes are concluded 
to be capable of changing the threshold value of transition from 
synchronous to asynchronous behavior among neurons [39]. 

Postnov et al. proposed a model containing three units (the 
presynaptic, and postsynaptic neurons and the glial cell) [5] . Their 
model can predict the long-term potentiation of the postsynaptic 
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Figure 1. Schematic of the three-unit model. N1: pyramidal cell; N2: interneuron; A: astrocyte. 
doi:10.1371/journal.pone.0080324.g001 



neuron. In the modeling study of astrocyte-neuron interaction, 
pyramidal cells and interneurons are often the focus [5,15,39,40]. 
Our study is likewise based the same coupled neurons. We will 
focus on the effect of astrocyte when information is transferred 
from pyramidal cell to interneurons, i.e., how the existence of 
astrocyte changes the response of interneuron to the firing pattern 
of pyramidal cell. 

Models 

Our model, which is schematized in Fig. 1, contains two 
conductance-based neurons and one astrocyte. Pyramidal cells are 
known to excite interneurons. By contrast, the interneurons inhibit 
the pyramidal cells. Thus, the two neurons in the model are 
coupled by excitatory and inhibitory synapses. The Hodgkin- 
Huxley equations have served a vital function in the theoretical 
understanding of neuronal behavior [41]. Following Ref. [26], we 
use the Hodgkin-Huxley equations to model the two neurons. The 
model equation describing the transmembrane potential contains 
sodium, potassium, and leak currents. The equations are given by 

C m = ~ gKn 4 J V x - Vx) -gm^hA V x - v Na ) 
-gdV x -v L ) + I asx + I ex + I sx , 

dm * n \ a 
—r- =%mA\ -m x )-p mx m x , 



= a«(l - n x ) - fi nx n x , 



—j-=ix ix (l -h x )-P hx h x , (1) 

where V x denotes the transmembrane potential of xth neuron 
(x= 1,2), and W represents the fraction of open Na + channels, and 
m x h x represents the fraction of open potassium channels. The 
values of parameters are listed in Table 1. The closing and 
opening rates of the gates are given by 

25 - V x -Vx 
°W = 0.1^5^ ,P mx = 4e is , 

e io -1 

an, = O.Ole-^r, p hx = — , 
e io +l 

10- V,. -Vx 
«« = 0.01 10 _ Vx ,^ = 0.125 e ^r, (2) 

e io + 1 

where I ex denotes the injected current input in xth neuron, Iasx is 
the feedback current received from the astrocyte by xth neuron, 
and I S x is the synaptic current received by xth neuron. Terman 
et al. [42] suggest that the neuron releases a neurotransmitter to 
the synaptic cleft depending on the membrane potential, and the 
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Table 1. Parameter values. 



parameter 


value 


c m 


1 uF/cm 2 


9k 


36.0 mS/cm 2 


Qua 


120.0 mS/cm 2 


9l 


0.3 mS/cm 2 


v K 


-12.0 mV 


V N a 


115 mV 


Vl 


1 0.6 mV 


0s 


85.0 


"s 


2.0 


<y- s 


0.1 


ft 


0.05 


g» 


0.1 


V S 1 


0.0 mV 


v se 


-85.0 mV 


Co 


2.0 \M 


Cl 


0.185 


Va 


6 s- 1 


Vb 


0.11 s -1 


Vc 


0.0 uM/s 


d, 


0.13 jiM 


d 2 


1 .049 |iM 


d 3 


0.9434 |*M 


d s 


0.08234 uM 


a 2 


0.2 uWT V 


(c 3 


0.1 flM 


Po 


160.0 nM 


TP 


0.00014 ms -1 


g ie 


variable 


0. 


variable 


The parameter values are obtained from references with slight modification. 
doi:1 0.1 371 /journal.pone.0080324.t001 



concentration of neurotransmitter released by xth neuron is given 
by 



T x 



1 



1 + e <.e s -v x )/<j s 



(3) 



Following Terman et al. [42] , the synaptic variable S x is introduced 
to explain the effect of the neurotransmitter release T y by jth 
neuron on the xth neuron, and the dynamic equation is given by 



Is2=gse(V 2 -v se )s l 



(5) 



where g se and g s , are the maximal conductance of the excitatory 
and inhibitory synapses. v se and v B - are the corresponding reversal 
potential. 

Astrocytes do not generate action potentials, i.e., the astrocytes 
are non-excitable electrically. The astrocytes respond to the 
neurotransmitter release in the synaptic cleft through IP3 
production [see Fig. 1]. Subsequently, elevation of IP3 concentra- 
tion induces the release of Ca 2+ from endoplasmic reticulum (ER), 
and then more Ca 2+ are released depending on the IP 3 -induced 
Ca 2+ elevation. The elevation of Ca 2+ above a certain threshold 
triggers the release of glial transmitters, which, in turn, will 
influence the dynamics of the neurons. We use Li-Rinzel model to 
describe the Ca 2+ exchange in the astrocyte [43] . This process 
contains three fluxes across the ER membrane: flux release 
through the ion channels (IP 3 Rs), removal of Ca 2+ by an ATP- 
dependent pump, and a leak. 



^= - cxVatnln^iC - C ER )- -jj^-^-civ h (C-C ER ) 



(6) 



dt 



--<x q (l-q)-P q q, 



with 



P C 
P + d x 



a q = a 2 d 2 , ; /L = a 2 C; 
n c 0 -C 

<~ER = 

Cl 



(7) 



where C denotes the Ca concentration in the intracellular 
space, q is the fraction of activated IP3R, and P is the IP 3 
concentration in the intracellular space. The values of parameters 
are listed in Table 1. The production of intracellular IP3 is 
modeled by 



dP 
dt 



Po-P 

TP 



(3) 



Nadkarni and Jung fit the experimental data [2 7] using the 



function of the current versus astrocytic Ca concentration 
I astro =2.1l®(lnc)lnc,c= C- 196.69(«M) 



(9) 



ds x 
dt 



tx s T y (l-s x )-[i s s x 



(4) 



Then, the synaptic current I sx received from each other by the 
two neurons in our model is 



Numerous physiological studies show that astrocytes release 
ATP, which has direct excitatory effects on hippocampal 
interneurons [44,45]. By contrast, astrocytes decrease pyramidal 
neuron excitability (Fig. 1) [36,46]. These findings suggest the 
following current I asx in Equation (1): 



Is\ =gsi(V\-v si )s 2 



Ias\ — — hi astro 1 1as 2 — + ^1 astro (10) 

where we introduce parameter acz(0,1) to account for the effect of 
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Figure 2. Time series of membrane potential in N1 and N2 for different parameter values. The successive spikes in N1 are induced by the 
injected current 7 e i = 10.0 ^A/cm 2 , and I e -,=0.0 /M/cm 2 by which the spikes can not be induced in N2; rp = 0.8 /jM/s. (a)A = 0, g„ = 0.5; (b)A = 0, 
gse = 0.9; (cU = 0.5, g„ = 0.5; (d)l = 0.5, g se = 0.9. 
doi:1 0.1 371 /journal.pone.0080324.g002 
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Figure 3. The critical value of g K for which the information is transferred from N1 to N2. I e2 = 0.0 /yA/cm 2 by which the spikes can not be 

induced in N2; j> = 0.8 /jM/s. 
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Figure 4. Time series of membrane potential in N1 and N2 for different parameter values. The values of parameters I e \, I t , 2 and r P are 

same as in Fig. 2. (a)(b)/. = 0.3, g. re = 0.9; (c)(d)/i = 0.5, g se = 0.9. Note that the time scales are much longer than that in Fig. 2. 
doi:1 0.1 371 /journal.pone.0080324.g004 



astrocytes. Given that X = 0, the effect of astrocyte is ignored, 
whereas when X — 1, the effect of astrocyte is considered fully. 

All parameter values are listed in Table 1 . We solve the model 
Equations (1) to (10) by using a fourth-order Runge-Kutta 
integration scheme with a time step 0.05, and simulations verify 
that further time step reduction does not significantly improve 
accuracy. 

Results and Discussion 

Ignoring the effects of astrocyte and synaptic current, i.e., 
lasx = 0 and I sx = 0, I ex > 6.24 fiA/ cm 2 is needed to generate 
persistent action potentials in the isolated H-H neuron. To study 
the information transmission from Nl to N2, we let I c \ = 10.0 jxAl 
cm 2 , and I e % — 0.0 //A/ cm 2 , for which the persistent action 
potentials are generated in Nl, but cannot be generated in N2 
on its own. While the effect of the astrocyte is ignored (X = 0), the 
persistent action potentials are found in N2 for large coupling 
strength g se , that means the information implied in the action 
potentials of Nl are transferred to N2. Comparing Figs. 2 (a) and 
(b), g se = 0.9 is sufficient for the information transmission, but 
g se — 0.5 is not. The results of calculation show that the critical 
value of gse for the information transmission is 0.56. When the 
effect of the astrocyte is considered, i.e., X > 0, the results are not 
significandy changed by the astrocyte [Figs. 2(c) and (d)]. As an 
example, in Fig. 3, the critical values are calculated for different 
values of X. Both for I e i = 10.0 and 20 /lA/ cm 2 , the critical value is 
independent of X but varies with different values of I e \. 



Bursting-like Spikes 

We now focus on a longer time scale. In Fig. 4, the value of g se is 
0.9, for which the persistent action potentials in Nl are successfully 
transferred to N2. When X = 03, the successive action potentials 
are transferred. Notably, bursting-like action potentials are found 
both in the Nl and N2 for X = 0.5. Bursting-like spikes (BLSs) are 
extensively found in experimental and modeling studies. Cressman 
Jr. et al. have studied the influence of sodium and potassium 
dynamics on neuronal behaviors using a single neuron model 
containing the effect of the glial cell [16,17]. They found the BLSs 
in some parameter regions, and the glial cell serves an important 
function in the appearance of BLSs. However, compared with our 
model, the glial cell in Ref. [16,17] modulates neuronal behaviour 
behavior by removing excess potassium from the extracellular 
space. Postnov et al. have found the postsynaptic neuron response 
to the presynaptic neuron by bursting-like firing in their modeling 
work regarding the effect of glial cell, but the presynaptic neuron 
fires successively [37]. This finding differs from our results because 
the BLSs always appear in Nl and N2 simultaneously. Theoret- 
ically, the resting membrane potential during the bursting spikes is 
attributed to the inhibitory effect of the astrocyte to Nl. In Fig. 5, 
the time series of the calcium concentration in the astrocyte and 
total current in Nl(/i =I as \ +hl +Isl) are depicted to correspond 
to Fig. 4. The calcium concentration is oscillating. When C is 
larger than 196.69 nM, the increase of the inhibitory current 
lasi (negative) will cause I\ to decrease to a low level. Otherwise, 
when C is larger than 196.69 nM, the inhibitory current vanishes, 
and I\ approaches 10 jUA/cm 2 . The red dashed lines in Fig. 5(b) 
and (d) represent 6.24 ^A/cm 2 . Obviously, while I\ decreases to a 
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Figure 5. Time series of calcium concentration and total current in N1 corresponding to that of membrane potential in Fig. 4. The 

values of parameters / eq , / e2 and r P are same as in Fig. 2. (a)(b) ).= 0.3, g se = 0.9; (c)(d) X= 0.5, g se = 0.9. The red dashed lines indicate the value 6.24 
//A/cm 2 . 
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value less than 6.24 /nA/cm , the Nl will possess the resting 
membrane potential. In Fig. 5(b), although I\ decreases owing to 
the larger C, I\ is always larger than 6.24 //A/ cm 2 . As a result, the 
successive firing of Nl will not be stopped. In Fig. 5 (d), I\ 
decreases to values less than 6.24 juA/ cm 2 periodically. When is 
less than 6.24 /iA/cm 2 , the Nl possesses the resting membrane 
potential, and BLSs are produced. Then, the BLSs are transferred 
to N2 through the excitatory synapse. 

The effect of astrocyte serves an important function in the 
production of the BLSs. As previously mentioned, the excitatory 
coupling strength determines the the information transmission 
from Nl to N2 significantly. Thus, we will identify the region of 
parameter X and g se , in which the BLSs are produced. 
Additionally, the IP 3 production rate rp has been proven to be 
associated with the expression level of mGluRs in astrocytes. The 
enhanced production of IP :i corresponds to over-expressed 
mGluRs. Over-expression of mGluRs has been reported to 
facilitate the seizure-like oscillations in the neurons [26]. In our 
study, three typical values of fp, 0.4, 0.5 and 0.8 /uM/s, are 
selected to represent the normal, intermediate, and enhanced 
expression level of mGluRs, respectively. The shadow regions in 
Fig. 6 are the parameter regions in which the BLSs can be found. 
First, the BLSs appear for an intermediate value of g se - Extremely 
large or small g se both make the calcium concentration approach a 
steady value less than 196.69 nM. In our model, the astrocyte fails 
to feedback to the neurons by I m iro when C is less than 
196.69 nM. Thus, the BLSs are not produced. Second, only if X 
is larger than a critical value do BLSs appear. Thus, we can 
conclude that the existence of astrocyte is an important condition 



for the production of the BLSs. Finally, the area of the shadow 
regions decreases sharply with decreasing rp. Enhanced expression 
level of mGluRs favors the BLSs. Although the models are 
different in previous literatures, the similar results have been 
obtained that the calcium dynamics in the astrocyte strongly affect 
the neural activity [5,15]. 

The rate of occurrence of the BLSs is then calculated. In Fig. 7, 
the rate / is approximately 0. 12 s~ 1 and is not very sensitive to the 
parameters, once the values of the parameters are within the 
shadow regions in Fig. 6. More accurately, / is maximum, and 
remains constant in the center of the shadow regions, f decreases 
when the parameter values change from the center to the edge of 
the regions. Furthermore, / increases with the enhancement of the 
expression level of mGluRs. Although Cressman Jr. et al. have not 
investigated the effect of astrocyte on the rate of the BLSs clearly, 
Ref. [16] shows that the rate increases with the enhancement of 
glial strength, and the rate is at the scale from 0.01 s~ l to 0.1 s . 
This finding is in accordance with our results qualitatively. 

Time Delay and Information Distortion 

Synaptic transmission is widely accepted to involve time delay 
attributed to the signal propagation time [47]. Theoretically, 
neuronal models with time delay have received considerable 
attention. Delay-induced coherent oscillation [48] is found in 
neuronal network as well as in other coupled systems. Delay- 
enhanced synchronization [23,49] may be relevant for neuronal 
networks to establish a concept of collective information processing 
in the presence of delayed information transmission. Our recent 
works find that delay cooperating with diversity can induce fruitful 
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Figure 6. Parameter region in which BLSs are produced. The values of parameters I,\ and are same as in Fig. 2. The value of rp equals to (a) 

0.4 fiM/s; (b) 0.5 //M/s; (c) 0.8 /jM/s. 

doi:10.1371/journal.pone.0080324.g006 



synchronization transitions [22]. Herein, the delay in the 
information transmission between two neurons will be verified in 
the presence of of astrocytes. As an example, in Fig. 8, the time 
series of V\ and Vj are recorded to show delay in the information 
transmission from Nl to N2. The spiking times in N2 always lag 
behind that in Nl. The time delay x is the time interval between 
two closest spikes in the two neurons. No matter whether the BLSs 
are produced or not, the time delay does exist in the information 
transmission. This time delay are also found in the previous 
modelling work studying the effect of astrocytes in neuron system 
[40] . Fig. 8 (b) and (d) show the time delay % corresponding to (a) 
and (c), respectively, x is not constant but oscillates irregularly. 
However, the oscillatory amplitude is not large; T possesses low- 
amplitude changing around a average value. 



The average value of x is calculated for different parameter 
values. Fig. 9 shows that with increasing g se , x decreases to a 
minimum first and then increases to a saturated value. The 
decrease of x for small g se corresponds to the increase of 
synchronization in Ref. [40]. The intermediate value of g se 
corresponding to the minimum x is about 2.96, and this value is 
independent of the expression level of mGluRs. Furthermore, for 
low expression level of mGluRs (rp = 0.2), x is totally independent 
on the value of X. With increasing rp, the minimal x will be 
influenced by X. Fig. 9(b) and (c) show that the minimal x reaches a 
maximum for an intermediate value of X, which is similar to the 
phenomenon of resonance found in random systems. 

Then, we will turn to another interesting phenomenon implied 
in Fig. 2 and 8. To exhibit this phenomenon clearly, the time series 




Figure 7. Rate of occurrence of the BLSs vs. the values of parameters g„ and X. The values of parameters I e \, / ( , 2: and r P are same as in Fig. 
6. 

doi:1 0.1 371 /journal.pone.0080324.g007 
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of V\ and V2 for different parameter values are depicted in Fig. 10. 
Notably, N2 does not respond to every spike in Nl through a 
corresponding spike accurately, i.e., large amounts of spikes are 



"missed" during the transmission from Nl to N2. Generally, the 
neuronal information is deemed to be coded in the spike timing or 
rate. Thus, the missing of spikes may relate to the distortion of the 
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Figure 9. Average rvs. the values of parameters g„ and X. The value of r P equals to (a) 0.2 /jM/s; (b) 0.5 /jM/s; (c) 0.8 fiM/s. 
doi:1 0.1 371 /journal.pone.0080324.g009 
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information transmission. Herein, we define the distortion ratio d r 
by the ratio between the spike number in Nl and N2. Obviously, 
all the spikes in N 1 respond by spiking in N2 for sufficient coupling 



strength. If g se is reduced, the number of missing spikes increases, 
i.e., d r increases. In Fig. 10, the values of the parameter X and Tp 
are set as 0.5 and 0.8, respectively. For these parameter values, the 
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BLSs are produced. In fact, the calculation shows that this kind of 
missing spike may occur whether the BLSs are produced or not. 

The distortion rate d r is calculated for different parameter 
values. Fig. 1 1 shows that the missing spike occurs mainly for small 
coupling strength g se . d r decreases to zero sharply if g se is 
increased from 0.56. The critical value of g se , for which d r 
decreases to zero, is about 1.06. Comparing the three figures in 
Fig. 11, d r is almost independent of rp and ).. Even the critical 
value oi gs e 1.06 does not change with the changing of rp and A. 
Thus, we can conclude that the effect of astrocyte does not serve 
an important function in the occurrence of missing spike. To 
exhibit the slight effect of astrocyte, the values of d r are amplified 
in Fig. 12. For intermediate or enhanced expression level of 
mGluRs, d r is non-zero with an intermediate value of )., whereas 
g se is larger than the critical value 1.06. We conclude that the 
effect of astrocyte induces the occurrence of miss of very few spike. 
Obviously, the accidental miss of spike does not result in the 
distortion of the information. 

Conclusions 

In this paper, the information transmission between neurons is 
studied by using a model that contains two neurons and one 
astrocyte. First, we identify the parameter region in which the 
information can be transferred from Nl to N2. The effect of 
astrocyte does not influence this parameter region. Secondly, in 



the parameter region for information transmission, we find BLSs 
in two neurons simultaneously. The parameter values for the 
occurrence of BLSs are also identified, and the results show that 
the higher expression level of mGluRs and the existence of 
astrocyte facilitate the occurrence of BLSs. Meanwhile, the rate for 
the occurrence of BLSs is calculated, and the rate is not very 
sensitive to the parameters. Third, time delay in information 
transmission is studied. The results show that T is not constant but 
oscillate with small amplitude. The average value of x is dependent 
on g se sensitively, but almost independent of rp and 1. Finally, we 
found amounts of spikes are "missed" during the transmission 
from Nl to N2. This distortion occurs mainly for small coupling 
strength g se . Although the astrocyte also induces very few missing 
spikes, it does not result in the distortion of the information. 

Although glial cells have been widely accepted to serve an 
important function in synaptic transmission in neuron system, 
theoretical knowledge on the mechanism of interaction between 
glial cell and neurons is lacking. The modelling studies in this 
paper can help us to understand the mechanism by which the 
astrocytes participate in neuronal information transmission. 
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